Nanotube Piezoelectricity 
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We combine ab initio, tight-binding methods and analytical theory to study piezoelectric effect 
of boron nitride nanotubes. We find that piezoelectricity of a heteropolar nanotube depends on 
its chirality and diameter and can be understood starting from the piezoelectric response of an 
isolated planar sheet, along with a structure specific mapping from the sheet onto the tube surface. 
We demonstrate that coupling between the uniaxial and shear deformation are only allowed in 
the nanotubes with lower chiral symmetry. Our study shows that piezoelectricity of nanotubes is 
fundamentally different from its counterpart in three dimensional (3D) bulk materials. 

PACS numbers: 77.65.-j,73.22.-f 



The physical properties of a nanotube along its ex- 
tended direction are controlled by the boundary condi- 
tions imposed along its wrapped direction. The exis- 
tence of both semiconducting and metallic forms of pure 
carbon nanotubes provides a striking example.^ The re- 
cently discovered electric polarization in heteropolar nan- 
otubes (e.g. Boron-Nitride) presents a new physical man- 
ifestation of this effect.^ Since the polarization can be 
modulated by elastic strains of the tube, these materi- 
als provide a new class of molecular piezoelectrics where 
mechanical strain is linearly coupled to an electric field. 
Piezoelectric nanotubes thus hold promise for application 
in nanometer scale sensors and actuators. 

In the modern quantum theory of polarized solids 
the electric polarization is computed from the geometric 
phase^ (Berry's phase) accumulated by the occupied elec- 
tronic states as one introduces a potential that adiabati- 
cally connects an unpolarized and polarized state of the 
system. For a BN nanotube the Berry's phase and hence 
the polarization is controlled by the periodic boundary 
condition on electronic wavefunctions.^ Piezoelectric ef- 
fect, on the other hand, is determined by the dependence 
of macroscopic polarization on the local strain induced 
effects: redistribution of the valence charge density, cur- 
vature induced rehybridization of the electronic orbitals 
and relaxations of the positions of the atoms on the tube 
walls which are all short range in character. Here we show 
that it is this latter character that allows the piezoelec- 
tric response to follow a simple transformation rule when 
the structure changes from a sheet to tube geometry. 

The prototypical example of piezoelectric nanotubes 
is found in the family of BN nanotubes where the al- 
ternation of group III(B) and group V(N) elements on 
the honeycomb lattice lowers the symmetry. A BN nan- 
otube can have a nonzero electric polarization^ unlike its 
planar counterpart where this is forbidden by the three- 
fold rotational symmetry of an isolated two dimensional 
(2D) BN sheet. However an elastic coplanar deformation 
of a BN sheet lowers its lattice symmetry, redistributes 
the valence charge and produces a nonzero polarization. 
Fig.Qillustrates the effect of a uniaxial strain {rjxx) and a 
shear strain {rjyx) of the BN sheet, with both distortions 
greatly exaggerated for clarity. These distortions induce 



FIG. 1: BN flat sheets under uniaxial strain r/xx (a) and shear 
strain rjxy (b). In both cases, threefold symmetry is broken 
and charge redistribution gives rise to a net dipole moment. 
The corresponding polarization directions (P) are marked by 



the electric dipole moments denoted by the arrows. The 
linear response of the electric polarization Pi to an ap- 
plied strain rjjk is described by the third rank piezoelec- 
tric tensor Cijk = dPi/dr]jk- The 3m symmetry of the 
unstrained sheet requires that the piezoelectric tensor is 
unchanged by threefold rotations of the lattice and the 
elements of the piezoelectric tensor obey the symmetry 

relation Cxxx — ^xyy — ^yxy — ^yyx' 

To study the microscopic origin of this behavior we 
first carried out ab initio calculations of the piezoelectric 
constants of the flat BN sheet using a planewave pseu- 
dopotential method based on density-functional theory 
within the local density approximation. The calculation 
is performed with the Abinit package^ using Troullier- 
Martins pseudopotentials^ with an energy cutoff of 45 
Hartree and 4x4x1 /c-point grid throughout. To create 
a computational cell that is periodic in all three spa- 
tial dimensions we stacked the BN sheets with an inter- 
planar distance 20 Bohr so that there is negligible wave 
function overlap between layers. The electronic polariza- 
tion was computed for a series of strained lattices using 



the Berry's phase formulation^ (discretized on a dense 
/c-point grid along the direction of the polarization) and 
the piezoelectric constants were obtained by calculating 
lattices with strains in the range of 1% < rjjk < 5%. 
By defining the positive direction to be the bond direc- 
tion from B atom to N atom as shown in Fig. ^ we find 
^xxx = — 0.12e/Bohr. Our calculations for the sheets 
with a shear strain rjxy have explicitly verified the symme- 
try relations among the piezoelectric tensor elements and 
showed that the piezoelectric properties of the sheet are 
controlled by a single coefficient. To carry out a system- 
atic study of the piezoelectric behavior of a large family 
of wrapped structures parameterized by integer indices 
(m,n)/ we combine the ab initio DFT method with a 
computationally less intensive albeit less accurate tight- 
binding (TB) method. We used a non-orthogonal basis 
set with four orbitals per site to describe the 2s and 2p 
atomic orbitals.^ We were able to benchmark our tight- 
binding method by comparing calculations of the piezo- 
electric constant of the BN sheet using both theories. 
We find that the TB theory yields e^xx = — 0.086 e/Bohr 
which is smaller than the ab initio result, though in ac- 
ceptable agreement. 

Note that the value of piezoelectric constant of a 2D 
BN sheet has a dimension charge per unit length. A 
quantitative comparison of the piezoelectric constant of 
the fiat sheet to the piezoelectric coefficients of 3D bulk 
material requires specification of the interlayer spacing 
and packing. For example, if we convert the above 2D 
fiat sheet value into a conventional "bulk" piezoelec- 
tric constant using the primitive interlayer separation of 
0.34nm,^ we find csd = 0.76 C/m^. This value is simi- 
lar in magnitude to 633 = 0.73 C/m^ of wultzite nitrides 
(e.g. GaN)^ and is larger than 0.12 C/m^ of piezoelectric 
polymer polyvinylidene fiuoride(PVDF-TFE).^^ Alterna- 
tively the piezoelectric constant computed from different 
dimensional system can be expressed as total dipole per 
stoichiometric unit. Using this convention, we find that 
the piezoelectric constant is 1.67 dipole/unit for the BN 
sheet, smaller than 1.98 dipole/unit for GaN wultzite. 
As we show below, when a sheet is wrapped to form a 
tube, there is a unique well defined relation between the 
piezoelectric constants of a tube and fiat sheet. 

Two high symmetry families of nanotubes are the zig- 
zag structures with wrapping indices (n, 0) and the arm- 
chair structures with wrapping indices (n,n). The one 
dimensional(lD) piezoelectric constants are defined as 
en = dPz/drjs and 614 = dPz/drjt^ where is the 
dipole moment per unit length and the (z, 5, t) indices 
in the tube frame refer, respectively, to the tube axis {z) 
and the uniaxial (5) and torsional {t) strains. In Fig. [3 
the top panel shows the structures of two representative 
small radius members of each family and in the bottom 
panel we plot their piezoelectric constants in unit com- 
parable with the two dimensional piezoelectric constants. 
Note that the ID piezoelectric constant is proportional 
to the tube circumference C through cid = C 620- We 
find that zigzag tubes exhibit a longitudinal piezoelectric 
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FIG. 2: Schematic structures of (a) stretched (n, 0) nanotubes 
and (b) twisted (n, n) nanotubes where arrows display the 
strain deformation direction in the tangential plane. Panels 
(c) and (d) show the calculated tube piezoelectric constant 
eii/C and 614/ C as functions of 1/n^ where C is the tube 
circumference. Contributions from a and tt electrons to the 
total piezoelectric response are separated. 

response for the case of uniaxial strain (extension or com- 
pression) but not for torsion. In contrast the armchair 
tubes have an electric dipole moment linearly coupled to 
torsion, but not to a uniaxial strain. The complementary 
strains, i.e. torsion for the zigzag structures and stretch 
for the armchair structures, produce a purely azimuthal 
dipole that integrates to zero on the surface of the cylin- 
der. 

For large radius tubes, one expects a correction to the 
piezoelectric constants of tube from its curvature, pro- 
portional to the inverse square of the tube radius. This 
can be seen in Fig. [21 bottom panel where we quantify 
this scaling behavior by plotting the calculated piezo- 
electric constants as a function of 1/n^. The data show 
that the tube piezoelectric constant rapidly approach the 
fiat sheet values with this scaling relation, but also that 
curvature effects remain quite small even for relatively 
small radius tubes. The data also show that the tt and 
a valence electrons have the same sign for both families 
of structures, with the tt electrons dominating the piezo- 
electric response, accounting for approximately 80% of 
the total. 

For a chiral tube the wrapping vector does not lie along 
a high symmetry direction of the 2D honeycomb lattice. 
This leads to a large and low symmetry translational unit 
cell for the chiral nanotube making a direct calculation 
of its piezoelectric properties cumbersome. We make use 
of the results for the high symmetry armchair and zigzag 
structures to develop an accurate scaling theory of the 
piezoelectric response of chiral tubes. Ignoring the finite 
radius corrections arising from the tube curvature the 
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FIG. 3: Piezoelectric response as function of the chiral an- 
gle in a sample of chiral nanotubes experiencing the uniaxial 
strain (a) and the shear strain (b). Solid black curves are 
the analytical result, Equation (0. Panel c shows ratio of 
the piezoelectric constants of chiral nanotubes to their flat 
sheet values plotted as a function of the inverse square ra- 
dius. The two branches are for the uniaxial en and torsional 
ei4 response. 



elements of the piezoelectric tensor are specified by ro- 
tating the known piezoelectric elements of the flat sheet 
onto the symmetry axes of the tube. Thus, defining the 
chiral angle as the angle between the axis of the tube 
and a 2D primitive translation vector we find 



^11 — C Cxxx 

ei4 



sin(3(9) 



C Cyxy cos{30). 



(1) 



Finite radius corrections to the predictions of Equation 
(PP) can then be obtained by comparing the results of this 
mapping to the values obtained from TB calculations on 
a selected set of chiral structures. Fig.|3Ji and b shows the 



result of this comparison, for uniaxial strain on the family 
of (5, m) tubes and for torsion on the families of (5,m), 
(6, m) and (12, m) families respectively. These data are 
very well described by the mapping of the 2D results, 
with surprisingly small corrections due to the tube cur- 
vature. The correction is quantified in Fig. |3t where we 
plot the ratio of the calculated piezoelectric modulus to 
its flat sheet (n oo) value as a function of for the 
family of (5, m) tubes. The deviations from the predic- 
tions of the flat sheet model are less than 15% over the 
entire range of structures we studied. 

Equation (QJ also reveals that the chiral tube has al- 
lowed linear piezoelectric coupling to both uniaxial strain 
and torsion, unlike the higher symmetry zigzag or arm- 
chair structures. Thus the long wavelength elastic en- 
ergy of a chiral tube generically has an anomalous cross 
term containing the product of the uniaxial and torsional 
strains. This implies that a tensile stress applied to a chi- 
ral tube induces torsion and conversely torsion induces a 
change of its length. Such a coupling is only possible for 
a chiral molecular structure, and indeed the coefficient 
of the cross term is a macroscopic manifestation of the 
underlying microscopic chirality of the nanotube. 

Recent progress in the synthesis of nanoscale ma- 
terials is demonstrating that many three dimensional 
lamellar phases can be fabricated in compact cylindrical 
structures. The appearance of pyroelectric and piezo- 
electric effects is a generic feature of these structures, and 
can be excluded only for special high symmetry wrap- 
pings. The methods we have developed and tested here 
for BN nanotubes should be widely applicable to study 
piezoelectric effects in this broader family of nanoscale 
materials. 
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